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ABSTRACT 

We present measurements of the optical/UV emission lines for a large homogeneous sample of 993 
quasars from the Large Bright Quasar Survey. Our largely automated technique accounts for continuum 
breaks and galactic reddening, and we perform multicomponent fits to emission line profiles, including 
the effects of blended iron emission, and of absorption lines both galactic and intrinsic. Here we describe 
the fitting algorithm and present the results of line fits to the LBQS sample, including upper limits to 
line equivalent widths when warranted. The distribution of measured line parameters, principally W\ 
and FWHM, are detailed for a variety of lines, including upper limits. We thus initiate a large-scale 
investigation of correlations between the high energy continuum and emission lines in quasars, to be ex- 
tended to complementary samples using similar techniques. High quality, reproducible measurements of 
emission lines for uniformly selected samples will advance our understanding of active galaxies, especially 
in a new era of large surveys selected by a variety of complementary methods. 

Subject headings: galaxies: active — quasars: emission lines — quasars: general — ultraviolet: galaxies 



1. INTRODUCTION 

Quasars, some of the most luminous objects in the Uni- 
verse, allow us to observe towards the start of cosmic time, 
probing extremes of physics, and illuminating the interven- 
ing matter. Though most of the QSOs now known were 
discovered by their unique spectral energy distributions 
(SEDs) or emission lines, several interrelated and funda- 
mental questions remain: 

1) What significant correlations exist between the ob- 
served SED and emission lines in large uniform quasar 
samples? 

2) How are the observed SED and emission lines affected 
by intrinsic absorption and orientation? 

3) Can the shape of the intrinsic SED be determined? 

4) What effect does the intrinsic SED have on emission 
lines and what docs that tell us about physical conditions 
in the QSO? 

The production of emission lines in QSO spectra is 
widely attributed to photoionization and heating of the 
emitting gas by the UV to X-ray continuum (e.g., Fer- 
land & Shields 1985; Krohk & Kallman 1988). Emission 
lines from a given ion are particularly sensitive to photons 
of energy above the corresponding ionization threshold, 
but ionization from excited states and heating via free- free 
and H^ absorption also determine the line's principal ion- 



izing/heating continuum. Many important lines respond 
to the extreme ultraviolet (EUV) or soft X-ray contin- 
uum. Unfortunately, the EUV band is severely obscured 
by Galactic absorption, although some constraints on the 
EUV ionizing continuum are available through analysis of 
the adjacent UV and soft X-ray windows. 

The overall similarity of QSO emission line spectra has 
on occasion been taken as evidence of fairly uniform phys- 
ical conditions in the broad emission line region (BELR). 
At first, this encouraged the assumption that clouds in the 
BELR inhabit a narrow swath of parameter space in den- 
sity, size, and ionization parameter ([/). Observed correla- 
tions of line equivalent widths (hereafter, Wx) with SEDs 
present challenges to geometric and photoionization mod- 
els that need to be answered. Photoionization pioneers 
such as Mushotzky & Ferland (1984) ran models on a sin- 
gle cloud, while photoionization models assuming a dis- 
tribution of clouds (i.e., not a single U) predict a strong 
dependence of W\ on the (assumed power law) slope of the 
SED (Dinette et al. 1989; Baldwin et al. 1995; Korista, 
Baldwin, & Ferland 1998). 

Large uniform samples of QSO emission line measure- 
ments are key to test photoionization models, but are also 
necessary for other important tests of quasar triggering 
and evolution. Tests of the binary vs. lens hypothesis 
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for quasar pairs at the same redshift require quantitative 
measures of the probability of finding the measured spec- 
tral similarity derived from complete samples (Mortlock, 
Webster & Francis 1999). Since optically-selected QSOs at 
similar redshift may tend to be more similar than 2 QSOs 
chosen entirely at random, studies of QSO spectral prop- 
erties should be based on samples large enough to provide 
statistically significant subsamples across a wide range of 
luminosity and redshift (Peng et al. 1999). 

In a small, complete sample of optically-selected QSOs 
studied by Laor et al. (1997), the strongest correlation 
found between X-ray continuum and optical emission line 
parameters was of the soft X-ray spectral slope {ax), and 
the FWHM of H/3. Strong correlations between ax and 
the relative strengths of [O III] A5007 and iron emission 
were seen both there and in previous studies, particularly 
Boroson & Green (1992) (hereafter BG92). A principal 
component analysis (PCA) was employed by BG92 to ex- 
amine an unwieldy set of emission line and continuum 
correlations within a sample of 87 QSOs taken from the 
Palomar-Green Bright Quasar Survey (BQS; Schmidt & 
Green 1983). This indicated that the strongest sources 
of variance in the optical spectra, the "primary eigenvec- 
tor" , involved an anti-correlation between the strength of 
[O III] A5007 and iron emission that was also linked to the 
FWHM of H/3. This work was extended by Marziani et 
al. (1996) (hereafter M96) and Sulentic et al. (2000) to 
include a study of the spectral region aroimd C IV A 1549. 
These papers identified Narrow-line type 1 Scyfert galax- 
ies (NLSl) and steep spectrum radio galaxies as occupy- 
ing opposite ends of the primary eigenvector. Other work 
utilizing the high spectral resolution of the Hubble Space 
Telescope Faint Object Spectrograph {HST FOS) showed 
a link between Al III A1859, Si III] A1892, and C III] A1909 
emission line strengths and the primary eigenvector but 
suffered from the small sample size and a focus on NLSl 
or the low redshift radio quiet objects in the BQS (Aoki 
& Yoshida 1999; Kuraszkiewicz ct'al. 2000). 

The samples of AGN used by most authors are almost 
entirely taken from BG92 and so do not represent a statis- 
tically complete sample and may not be truly representa- 
tive of the QSO population (Wisotzki et al. 2000). They 
also deal with the strongest emission features in the OUV 
spectrum and, while several hypotheses have been pro- 
posed, the physical origin of the observed correlations re- 
mains a puzzle. There is also a concern that biases within 
these small samples may have exaggerated the significance 
of the primary eigenvector. 

Even in the clear presence of extrinsic effects such as ab- 
sorption along the line of sight, emission lines have been 
used to infer the strength and shape of the high energy 
SED. As an example, the similarity of UV emission-line 
properties in broad absorption line (BAL) and non-BAL 
QSOs (Weymann et al. 1991) has been cited as evidence 
that orientation is the cause of the BAL phenomenon (i.e., 
that all radio-quiet QSOs have BAL clouds). However, 
BALQSOs are now known to exhibit markedly weak X- 
ray emission as a class (Green & Mathur 1996; Gallagher 
et al. 1999; Brandt, Laor, & Wills 2000). But if similar 
emission lines indeed vouch for similar intrinsic high en- 
ergy SEDs, then the large X-ray to optical continuum fiux 
ratio {aox) values observed for BALQSOs are likely to be 



caused by strong absorption along the line-of-sight rather 
than by differences in their intrinsic SEDs. In fact, the 
soft X-ray deficit requires absorption at least an order of 
magnitude larger than those estimated from optical/UV 
(OUV) spectra alone and the UV and X-ray absorbers 
have yet to be definitively identified as one (e.g., see discus- 
sion in Mathur, Wilkes, & Elvis 1998). Since BALQSOs 
may be heavily absorbed, the question of whether simi- 
lar emission lines are testimony for similar intrinsic SEDs 
must be answered through the systematic study of line and 
continuum correlations in unabsorbed QSOs. 

Even in supposedly unabsorbed QSOs, a strong inverse 
correlation exists between optical iron emission, and O III 
A 5007 in QSOs (BG92). Iron emission is under-predicted 
by a large factor in standard photoionization models rela- 
tive to observations (Collin & Joly 2000), and the contin- 
uum source responsible for creating rest-frame optical iron 
emission lines comes almost entirely from X-rays above 
about IkeV. QSOs with large aox of which BALQSOs 
represent an extreme, have not only weak O III, but feeble 
narrow line emission in general (Green 1996). We suggest 
that the observed correlations may be an effect of absorp- 
tion of the intrinsic high energy continuum by thick clouds 
interior to the NELR. The X-ray heating of these clouds 
makes them emit copious Fe II, lines, but prevents the 
ionizing continuum from reaching the NELR. This is re- 
inforced by another observed correlation - that the most 
highly absorbed BALQSOs, the low-ionization BALQSOs, 
albeit X-ray quiet, are generally very strong optical Fe II, 
emitters (Stocke et al. 1992; Lipari, Terlevich, & Mac- 
chetto 1993). Our interpretation needs to be tested across 
a range of absorption, in large representative QSO sam- 
ples. 

We have therefore undertaken a major study of quasar 
line emission, with careful accounting for absorption lines 
and blended iron emission, using largely automated proce- 
dures on carefully-selected samples. The analysis of sam- 
ples of QSOs with more than a few himdred spectra re- 
quires some amount of automation to give consistent re- 
sults. The largest sample of QSO spectra currently avail- 
able is that of the Large Bright Quasar Survey and here we 
describe the initial results from the measurement of this 
sample and the analysis methods that will also be applied 
to other large samples of QSO spectra that will become 
available in the near future. 

2. THE LBQS SAMPLE 

The Large Bright Quasar Survey (LBQS; see Hewett, 
Foltz & Chaffee 1995) is a sample of 1058 QSOs selected 
using the Automatic Plate Measuring Machine from UK 
Schmidt photographic direct and objective prism plates. 
The combination of quantifiable selection techniques, in- 
cluding overall spectral shape, strong emission lines, and 
redshifted absorption features has been shown to be highly 
efficient at finding QSOs with 0.2 < z < 3.3, a signifi- 
cantly broader range than previous work. The LBQS thus 
avoids selection effects common in other optical quasar 
samples that tend to exchide weak-lined quasars or to 
under-sample certain redshift ranges or colors. There ap- 
pears to be a deficiency in the numbers of QSOs with 
z ~ 0.8 in the LBQS sample, although the effect is not sta- 
tistically significant (Hewett et al. 1995). The LBQS has 
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been for the last decade the principal source of interme- 
diate brightness optically-selected quasars available; 1/9 
of the Veron-Cetty & Veron (1998) catalog of QSOs and 
~ 50% of all known quasars within a similar magnitude 
range. Radio data are currently available for about 1/3 
of the sample (Hooper et al. 1995) and soft X-ray fluxes 
and upper limits for at least 85% of the sample can be 
obtained from the ROSAT AW-Sky Survey. Follow-up op- 
tical spectra with S/N«10 in the continuum at 4500A and 
6 — lOA resolution were obtained between 1986 and 1990 
at the MMT on Mt Hopkins, Arizona and the 2.5m duPont 
in Las Campanas, Chile for all LBQS QSOs. Error spectra 
(Icr) are available for 1009 of the objects. The selection cri- 
teria and observing procedures for each of the survey fields 
can be found in Hewett et al. (1995) and references therein 
(LBQS Papers I- VI). Three QSOs appear here that were 
identified subsequently to the 1055 QSOs in the Hewett et 
al. (1995) LBQS catalog. They are 0052-^0148 (z=0.595 
), 1027-0149 (z=0.754), and 2132-4227 (z=0.569), see 
Hewett, Foltz, & Chaffee (2001; in preparation) for fur- 
ther details. 

3. ANALYSIS OF THE SPECTRA 

The LBQS spectra were analyzed using Sherpa ^, a gen- 
eralized fitting engine designed primarily for spatially re- 
solved spectroscopy of observations with NASA's Chandra 
X-ray Observatory. Sherpa enables data to be modeled 
using a variety of optimization methods and with a num- 
ber of built-in statistical tests. The ease with which user- 
models can be created and the ability to simultaneously 
model a number of input data sets that may be in different 
formats (e.g. modeling an ASCII and a PHA style FITS 
file together) makes Sherpa suitable for our goals. In brief, 
the process begins with modeling the continuum emission 
and then accounting for any emission by iron complexes in 
the spectra. We then model the emission lines using Gaus- 
sian profiles and, after a search for significant absorption 
features, the results are inspected and the emission line 
fitting procedure repeated to improve the model of each 
spectrum. The model parameters were determined from a 
minimization of the statistic with modified calculations 
of uncertainties in each bin (Gehrels 1986). We found that 
a Powell optimization method gave a balance between an 
efficient processing time and consistency of results. Be- 
cause of the large database of observations the modeling 
of the spectra proceeds in a largely non-interactive man- 
ner. 

We do not present measurements of emission lines for 
objects with strong BALs. The diversity in the strength 
and shape of BALs precludes even semi-automated mea- 
surements, and their emission line properties have already 
been well-characterized for a strongly overlapping sample 
in Weymann et al. (1991). We also exclude objects that 
have multiple narrow absorption features that significantly 
mask the true shape and strength of the emission features 
of Lyman a and C IV A1549. A list of the 65 objects 
excluded from this study can be found in Table 1. This 
list includes two objects not previously classified as BAL 
in LBQS papers I-VI, 0059-2545 and 1242-f 1737. The 
measurements for the remaining 993 LBQS spectra are 
presented here. 



3.1. The continuum 

Our first step is to model the continuum emission by 
fitting one or two power laws to ranges in the spectrum 
that are free from any strong emission lines. We found the 
second power law component to be necessary due to the 
steepening of the continuum towards the UV in low red- 
shift objects, i.e. for spectra that extended redward of rest 
frame 4200 A. A list of the continuum modeling 'windows', 
along with the nearest strong emission lines, can be found 
in Table 2. We apply a Galactic reddening correction to 
the power law continua that follows the prescription given 
by CardeUi, Clayton, & Mathis (1989) which uses: 



A{X) = E{B -V)[aRY + b] 

where a and b are polynomial functions of wavelength de- 
rived for 3.3/im > A > 1000 A. We use a ratio of total to 
selective extinction of i?v = 3.1, the standard value for the 
diffuse ISM, and take the relationship between the color 
excess and the line of sight column of neutral Hydrogen 
for each object to be 



E{B -V) = 



A^H(Galactic)(102 
58l) 



(Bohlin, Savage, & Drake 1978). 

Many of the spectra in the LBQS sample depart from 
a simple reddened power law continuum for wavelengths 
blueward of ^ 4000A in the observed frame as noted in 
the original LBQS papers. This depression in the observed 
ultraviolet is due to the use of 2.5 arcsecond circular aper- 
tures in the MMT observations coupled with the effects 
of atmospheric dispersion at large zenith distances and a 
guider system that was red-sensitive. To model this part 
of the continuum, we use a polynomial of up to 2nd or- 
der with the flux tied to that of the power law continuum 
model in the nearest continuum window with Ac (observed 
frame) > 4000A. We did not apply a polynomial contin- 
uum to objects with z > 2.42 because the Lyman a emis- 
sion line appears redward of 4000A and the continuum on 
the blue side of Lyman a can be significantly reduced by 
Lyman a forest lines. 

The polynomial continuum model is not physically 
meaningful and may be inaccurate, particularly where 
strong iron emission complexes appear around the Mg II 
A2800 emission line, i.e. for 

2200A ^ Arcst <, 3100A 

We investigated the effect of the polynomial continuum 
on W\ measurements and on the error estimates for the 
emission line parameters by comparing lines of the same 
species that occurred in the polynomial continuum re- 
gion to those that did not. A Kolmogorov-Smirnov test 
(e.g.. Press et al. 1992) shows that the distributions of 
W\ measurements for the UV iron emission. Si IV -fO IV 
A1400, C IV, and Mg II emission lines are significantly 
different when measured above the polynomial continuum 
compared to when measured using a power law contin- 
uum. The median increase in W\ for C IV and Mg II 



^Sherpa is available from the Chandra science center { ittp : //asc .harvard. edu ) 
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is ^ 20% when these hncs fall above the polynomial con- 
tinuum but genuine observational trends like the Baldwin 
Effect (Baldwin 1977; Osmer, Porter, & Green 1994) may 
account for this measurement trend. The intrinsic con- 
tinuum luminosity of objects with C IV redshifted to A 
(observed frame) > 4000A is on average twice that of the 
lower redshift objects where C IV falls above the polyno- 
mial continuum. A 20% increase in Wx for C IV and Mg II 
is predicted using the relationships between Wx and con- 
tinuum luminosity determined by Zamorani ct al. (1992) 
but may not fully account for the differences seen in the 
measurements of Wx for these and other emission lines. 
Rather than correcting for this effect, we added a system- 
atic uncertainty of 20% to the estimation of emission line 
strengths and upper limits for all emission lines measured 
above the polynomial continuum. 

Finally, we visually inspected the continuum model for 
each observation, and found that minor adjustments were 
required in ~ 20% of the sample; a large proportion of 
these were caused by strong intrinsic absorption lines. 

3.2. Iron emission 

Our second step accounts for iron emission line com- 
plexes in the spectra. This emission occurs most strongly 
in the regions around the Mg II A2800 and H/3 emission 
lines. We follow the prescription used by BG92 of sub- 
tracting a template of iron emission lines created from the 
spectrum of I Zw 1, an NLSl which exhibits the typically 
strong iron emission of this class of AGN. We use two tem- 
plates: the optical emission line template is identical to 
that used by BG92 and corrects for iron emission between 
4400A < Arost < 7000A, while the UV template was devel- 
oped by Vestergaard & Wilkes (2000) from HST FOS ob- 
servations of I Zw 1 and covers the region 1250 A < Arest < 
3100A. The FWHM of the iron emission features seen in 
I Zw 1 is similar to that of the H/3 emission line (900 
km ) and so sets the minimum template FWHM. We 
convolved this template with a series of Gaussian func- 
tions of increasing width to produce a grid of 38 template 
spectra with 900 <; FWHM ^ 10000 km s~^ , conserving 
the total flux in each template. This provides us with a 
nominal resolution of 250 km s^^, more than adequate for 
the quality of the spectra in the LBQS sample. 

The first part in modeling the iron emission is to com- 
pare the 2000 km s"-'^ templates to windows in each spec- 
trum on both sides of the Mg II A2800 emission line and 
redward of C III] A1909 for the UV emission, and to ei- 
ther side of the H/3 + [O III] AA4959,5007 complex for the 
optical emission (see Table 2). After the relative ampli- 
tude of the iron emission has been modeled, the full grid 
of templates is applied to the spectra to determine the ap- 
proximate FWHM of the iron emission. This is followed by 
a final modehng of the template amplitude simultaneously 
with the FWHM of the template spectra. Adjustments 
are made from a visual inspection of the results to give the 
best representation of the iron emission complexes present 
in each spectrum. 

Note that we measure the strengths of the iron emission 
in the UV and optical independently. This is important 
as the relative strength of the UV and optical iron emis- 
sion may allow an approximate determination of the con- 
ditions in the emitting region, e.g. density and U (Verner 



ct al. 1999). UV and optical iron emission may also cor- 
relate differently with X-ray emission. Green et al. (1995) 
found that QSOs in the LBQS with strong UV FeH emis- 
sion (based on the iron feature under [Ne IV]A2423) are 
anomalously X-ray bright in the ROSAT passband. By 
contrast, several studies (e.g., Corbin 1993, Lawrence et 
al. 1997) find an anti-correlation between soft X-ray lumi- 
nosity and R{Fe II) (the equivalent width ratio of optical 
iron emission to H/3). Several such intriguing correlations 
have been noted, but among samples of varying size and 
using heterogeneous measurement techniques. Our project 
attempts to remedy this situation using large samples and 
uniform analysis. 

When first scaling the UV iron template in the LBQS 
spectra to the iron emission complexes straddling Mg II 
(~ AA2200 — 3300), we found that the template tended 
to over-predict the iron emission line strength blueward 
of C IV A1549. This occurred in 27 of the 238 spectra 
that cover both C IV and Mg II and that have measurable 
iron emission. The UV iron template is based on an HST 
FOS spectrum of I Zw 1 for which the region blueward of 
C IV was observed six months earlier than the remaining 
spectrum. I Zw 1 clearly brightened between the two ob- 
servations, for which reason the bluest sub-spectrum was 
scaled to match the redder spectrum. The reasons for the 
over-prediction of iron flux blueward of C IV in our sam- 
ple when using the I Zw 1 template are not clear, but may 
include a change in iron emission strength and multiplet 
ratios with luminosity indicating variations in the physi- 
cal conditions to which the iron emission is sensitive (e.g. 
Netzer 1980). See Vestergaard & Wilkes (2000) for further 
details and discussion. To compensate for this effect, the 
UV iron template flux in the region Aicst < 1530 A was 
reduced by 50%, resulting in an improved agreement with 
all the sample spectra. 

3.3. Emission lines 

The emission lines in the spectra are modeled by the ad- 
dition of Gaussian features at the expected wavelengths to 
the continuum -|- iron template model. The initial value of 
the peak flux near the position of the line is estimated from 
the raw data assuming a continuum level under each emis- 
sion line calculated by fitting a straight line between the 
continuum windows nearest the line to be modeled. The 
initial estimate of the FWHM of the emission line compo- 
nent is assumed to be 3000 to 5000 km s^^ for broad com- 
ponents, and 350 km s~^ for narrow components. These 
reflect the typical widths of emission lines produced in the 
BELR and NELR (e.g. Peterson 1997). Where the S/N 
in the spectra are high enough, the Lyman a, C IV, Mg II, 
and H/5 emission lines may be modeled using two Gaus- 
sian components. The initial value of the peak amplitude 
of each component is set to be 40% of the peak flux es- 
timated from the data. Lower S/N spectra are modeled 
with a single Gaussian but a visual inspection of the re- 
sults is made and a second Gaussian component is added 
where significant residuals appear. 

The FWHM and peak amplitude of the Gaussian model 
components are then optimized within a spectral region 
covering the emission line. The position of the emission 
line is fixed at the expected wavelength in the initial mod- 
eling until after the presence of narrow absorption features 
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has been determined (§3.4). The region of the spectrum 
used in modehng an emission feature is chosen to best op- 
timize the component of interest, e.g. the Lyman a emis- 
sion hne is simuhaneously fitted with N V A1240 and O I 
A1305 because the broad component to Lyman a may have 
a significant effect on the spectrum near OLA Ust of the 
emission hnes modeled in the LBQS spectra is presented 
in Table 3. There are a number of important emission 
lines (e.g. C II A1336) that are not modeled here due to 
the quality of the spectra, but will be included in the ap- 
plication of this procedure to samples of higher spectral 
resolution (e.g. HST FOS spectra). 

3.4. Absorption lines 

One of the goals of this study is to determine the fre- 
quency of occurrence of absorption in the spectra of QSOs. 
We used the program FINDSL (Aldcroft 1993) to detect 
significant narrow absorption features and model them 
with multi-component Gaussian fitting. Our input to 
FINDSL is the sum of the reddened continuum, iron tem- 
plate (s) and all the emission line profiles from the first 
Sherpa modeling. This constitutes a 'continuum' from 
which FINDSL detects significant deviations in the ob- 
served spectrum. We rescale the error array of each spec- 
trum for use with FINDSL so that xt becomes unity, which 
yields a better estimate of the true error array and aids in 
the detection of significant absorption features. We ex- 
clude the Lyman a forest region blueward of rest frame 
1065A from the absorption line detection, and also ex- 
clude the Balmer continuum region (between 3360A - 
3960A ) where the global continuum model may lie below 
the spectrum causing many spurious absorption lines to 
be detected. 

The absorption line measurements from FINDSL were 
then included with the emission line and continuum re- 
sults for further modeling using Sherpa. In this second it- 
eration the emission line positions are also fitted, and the 
absorption lines modeled simultaneously with the emis- 
sion line components. The results from this automated 
process were inspected and adjustments made to those 
spectra that the algorithms could not successfully model. 
This was not unexpected as the wide variety of spectra 
in the sample and the chance superposition of absorption 
features on the emission lines makes it unlikely that any 
fully automatic procedure would be 100% successful. 

3.5. Error analysis 

An estimate of the 2a error range for each parameter 
of the continuum and emission line components was de- 
termined from the confidence interval bounds (A^^ = 
4.0). Where the amplitude of an emission line could not 
be constrained to within 2(7, the line position was reset 
to the expected wavelength and the FWHM to a median 
value that was determined from the distribution of well 
constrained emission line measurements within the LBQS 
sample. These fixed emission line parameters are given in 
Table 3. The emission line amplitude was then re-fitted 
and, if remaining unconstrained at the 2a level, then the 
2a upper limit on the amplitude was estimated. The de- 
termination of upper limits for lines that may be present 
at low flux levels in the spectra of QSOs is one of the pri- 



mary objectives of this work. This allows the use of sur- 
vival analysis techniques in determining a more realistic 
distribution of emission line strengths. 

4. EMISSION LINE RESULTS 

The total number of emission lines measured from 
the spectra in the LBQS sample are presented in Ta- 
ble 4. The numbers of upper limits measured and the 
number of two component emission lines are also tab- 
ulated. The full tabic of all emission line measure- 
ments and uncertainties of the 993 spectra modeled here 
is available from the ApJ in the electroni c journal at 
http : // www ■ j ournals . uchicago . edu/ Ap J/ , but due to 
the large size cannot be reproduced in the printed jour- 
nal. To show the form and content of the large electronic 
table we present in Table 5 a digested version containing 
measurements of the emission lines present in three LBQS 
spectra. The emission line measurements are quoted to 
a significance level determined by the 2a uncertainties in 
the parameters (e.g. Bevington & Robinson 1992). The 
format is as follows: Each object is represented by a row 
containing the designation and redshift followed by a row 
for each measured emission line. The name of the emission 
line is given in column (1), column (2) gives the FWHM in 
km s~^ , and column (3) gives the offset of the peak of the 
Gaussian emission line model, in km s^^ , from the ex- 
pected position based on the tabulated redshift. Note that 
no peak offset measurements appear for the iron emission 
line measurements. Column (4) gives the rest frame equiv- 
alent width of the emission line in A. Each emission line 
parameter is quoted with positive and negative 2a error 
estimates. The errors quoted for W\ are based in the un- 
certainties in the amplitude and FWHM of the Gaussian 
model and do not include an error from an uncertainty 
in the underlying continuum flux level. For emission lines 
where only an upper limit on W\ could be measured, 
there are no values quoted for the peak offset because we 
flxed the position of the line at its expected wavelength. 
The FWHM value also has no associated errors in such a 
case because it was fixed at an approximate median value 
for the sample. (See §3.5.) Note that some detected but 
poorly constrained lines may have an error estimate for 
the W\ measurement even while the FWHM is flxed (see 
§3.5 and Table 4). Finally, column (5) gives the number 
of narrow absorption features used in modeling the emis- 
sion lines. They are tabulated in the row of the closest 
emission line to the position of the absorption feature, e.g. 
see Lyman a (broad) and O I for LBQS 2354-0134 in Ta- 
ble 5. The format of the electronic version of Table 5 is 
explained in the Appendix. 

To generate a measure of the strength of iron emission in 
each spectrum, the integrated flux of the model template 
in the regions 2240A - 2655A (blueward of Mg II) or 
4434A - 4684A (blueward of H/3) is used and combined 
with the continuum specific flux at 2448A and 4559A re- 
spectively to calculate W\. For objects where these regions 
were not present in the observed spectrum, we extrapolate 
the continuum to these positions. Note that we integrate 
the optical iron template flux over a wavelength range that 
matches the window used by BG92. 

We present in Table 6 the continuum parameters for the 
same LBQS spectra that appear in Table 5. The full table 
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is available in electronic form with one row for each QSO in 
the sample. The format of Table 6 is as follows. Columns 
(1), (2) and (3) give the QSO designation, redshift, and 
the Galactic neutral Hydrogen Column density (in units 
of 10^° cm~^) respectively. The latter measurements are 
taken from the Bell Laboratory H I survey (Stark et al. 
1992). Columns (4) - (6) contain the measured values 
of the power law continuum; continuum slope, the rest 
frame wavelength at which the continuum is normalized, 
and the continuum flux at that wavelength (ergs cm~^ 
A~^) respectively, along with their 2a uncertainties. Note 
that in Table 6 the normalization units are 10~^^ but in 
the electronic table the units are 10^"'^^. Column (7) lists 
the slope of the second power law continuum if used, the 
normalization of which is identical to the first power law 
continuum. Column (8) presents the order of the polyno- 
mial continuum if required to model the blue region of the 
spectra. The format of the electronic version of Table 6 is 
explained in the Appendix. 

We present in Figures 3a, 3b, and 3c the profiles and 
residuals from the models of the three spectra presented 
in Tables 5 and 6. The measurements presented here were 
chosen to show a range of spectral qualities present in the 
LBQS sample. Fig. 3 includes a number of panels for 
each spectrum in each case all wavelengths are observed 
frame and fltixes are 10"^'' ergs cm~^ s~^ A~^. The top 
panel has three sections, the upper section shows the con- 
tinuum model overlying the spectrum, including error bars 
on each bin, the middle section displays the iron emission 
template plotted over the spectrum (no error bars) and 
the lower section shows template profile alone. Note that 
the flux scales arc different in these sections to display the 
template emission more clearly. 

This panel is followed by a number of separate pan- 
els, one for each emission line region. These emission line 
panels have three sections; the top section shows the to- 
tal continuum + emission line model over the spectrum. 
(The model will include flux from the iron emission line 
template(s) if used. For clarity, flux errors on each bin are 
not shown in this section.) The middle section shows the 
residuals from the emission line model and include the er- 
ror bars for each spectral bin. The lower section shows the 
profiles of the emission line models, with multiple compo- 
nents separated into single Gaussian profiles. The shape 
of the iron emission template used in the modeling of the 
spectrum is also shown in the lower panel. Dashed vertical 
lines indicate the expected position of emission lines based 
on tabulated redshifts. Absorption line profiles are shown 
along the top edge of the lower panels, e.g. the Lyman a 
region panel for LBQS 2354—0134 in Fig. 3c. For the pan- 
els that show the emission line models of [Nc V] , [O II] , and 
[Ne 111] the top sections also plot the local continua used 
for these lines (see Appendix). Note that to help visual- 
ize the quality of the emission line models the wavelength 
scale is identical for each of the emission line panels and 
within each panel the flux scale is identical in the sections 
displaying the total model, residuals, and individual model 
components. 

5. OVERALL STATISTICS 

The analysis of the LBQS sample of QSO spectra has 
provided us with over 8000 emission line measurements, of 



which approximately 30% are upper limits to the equiv- 
alent width of low intensity lines. Without the correct 
inclusion of these "censored" data in the analysis of this 
complete sample, a biased estimation of the properties may 
occur and a number of important relationships may be 
masked or, even worse, appear more significant than in 
reality. 

We use a non-parametric survival analysis technique 
to estimate the means and medians that characterize the 
emission line parameter distributions. The Kaplan-Meier 
(KM) estimator of a sample distribution is a maximum 
likelihood method that provides a reconstruction of infor- 
mation lost by censoring of data. This well established sta- 
tistical technique has been throughly examined for use in 
astronomical studies by Feigelson & Nelson (1985), Isobe, 
Feigelson, & Nelson (1986) and references therein. We 
recommend these papers to the reader. 

The astronomical survival statistics package ASURV 
Rev 1.1 (Isobe & Feigelson 1990; LaValley, Isobe & Feigel- 
son 1992) was used to provide the means, errors on the 
mean, and the medians of the Wx distributions for the 
LBQS sample given in Table 7 and the distributions of 
Wx shown in Figure 4. The format of Table 7 is as follows. 
Column (1) lists the emission line or line blend measured. 
The distribution for single Gaussian component models 
are tabulated separately from narrow and broad compo- 
nents. The distribution of sum of the Wx of the broad 
and narrow component models included with the single 
component measurements is also tabulated. Columns (2) 
to (4) give the number of detected emission lines, their 
mean Wx and the standard deviation (SD) of the distri- 
bution of Wx ■ Columns (5) and (6) list the total number 
of emission lines measured for each species and the num- 
ber of upper limits estimated for Wx- Columns (7) and 
(8) give the KM means, error on the means, and medians 
of Wx for each line. Columns (9) to (11) give the number, 
mean and median of the FWHM of the Gaussian compo- 
nents used to model each emission feature. All equivalent 
widths are rest frame and the error on the Wx means is 
estimated from the Icr of the KM reconstructed distribu- 
tion. Note that there are no upper limit measurements 
for two component emission line fits, in cases where any of 
the parameters of the two components could not be con- 
strained to 2a, a single Gaussian was used to model the 
line. 

We should add a note of caution here about the recon- 
structed distribution of emission line properties for sam- 
ples that contain a large proportion of censored data (e.g. 
[Ne V], He II A4686). Formally, an emission line is de- 
tected if the amplitude of a feature can be constrained 
to within 2a, based on the confidence interval boimds. 
The performance of the KM estimator degrades if the cen- 
sored fraction is high or if the pattern of censored mea- 
surements is not random. We do not believe the latter is a 
concern here because the observing procedure used in the 
generation of the sample spectra produced spectra of sim- 
ilar continuum S/N, independent of the strength of emis- 
sion features and the intrinsic luminosity of the QSO. In 
support of this conclusion consider the statistics of the O I 
emission line, where approximately half of the spectra with 
coverage in the region of this line possess detected emission 
features. Detected emission lines have 0.6 ^ Wx ^ 14.0 
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A but the range of upper limits for the undetected hnes 
is 0.2 < Wx < 15. OA . The large proportion of censored 
data for some emission lines lowers our confidence in the 
KM estimate of the means and medians of a reconstructed 
W\ distribution. Any conclusions based on these emission 
lines should rather use the results from the detected lines 
alone. 

The histograms shown in Figure 4 are taken from the 
KM estimation of the number of data points in each bin. 
The W\ (A) and FWHM ( km s~i ) distribution for each 
emission line are shown in the upper and lower panels for 
each emission line respectively. For emission lines that 
have been modeled using two Gaussian components, we 
include the distribution for the narrow and the broad com- 
ponents separately from the single component lines. For 
completeness we also show in Figure 4 the distribution 
of the full sample of measured lines where the W\ of 
broad+narrow components of an emission line have been 
summed and included in the single component distribution 
(dashed line histograms). 

6. DISCUSSION 

To examine the robustness of the analysis procedure in 
OUV spectra with improved S/N and spectral resolution 
we have modeled two composite QSO spectra using the 
techniques described above. The first was constructed by 
Francis ct al. 1991 (hereafter F91) from 718 spectra in 
the LBQS sample itself, and the second from HST FOS 
observations of 101 QSOs presented in Zheng et al. (1997) 
(hereafter Z97). We find that even where the limit of two 
Gaussian components for an emission line does not repro- 
duce the line profile very well, the measured Wx is within 
10% of values tabulated in the above papers. The excep- 
tion is for Mg 11 where the F91 measurement is higher and 
the Z97 measurement is lower than that measured using 
the techniques presented above. Both discrepancies can be 
attributed to the presence of iron emission features blend- 
ing into the wings of the Mg II emission line. Although F91 
fitted a global continuum and corrected for iron emission 
using local spline fits to emission complexes, the Wx was 
measured as the total flux above the continuum and iron 
emission across more than 50A in the spectrum and so 
may include some hidden iron emission. This is also likely 
to be the explanation of a similar discrepancy in the mea- 
surements of H7 + [O III] A4363 and [O III] A4959 in 
the F91 composite. The measurements in Z97 use a lo- 
cal continuum which is difficult to determine accurately 
in regions where iron emission complexes are strong, and 
may result in their lower Wx measurement. This confirms 
that correct accounting for iron emission is crucial in the 
measurement of emission lines in QSO spectra. 

We can also compare the global properties of the LBQS 
emission line measurements with the measurements of 
these composite QSO spectra. The Z97 composite shows 
larger Wx in Lyman /? + O VI, Lyman a, and C IV. 
The discrepancy remains even if the means of only our de- 
tected emission lines are compared. (The inclusion of the 
censored data will naturally lower the means of a distri- 
bution, but has little effect on the distribution of strong 
lines like Lyman a and C IV which are well constrained 
in nearly all cases). We suspect that the difference is at- 
tributable to the large (60%) fraction of radio loud ob- 



jects in the heterogeneous sample used to construct the 
Z97 composite. The LBQS sample is predominantly Ra- 
dio Quiet (RQ) with ~ 10% Radio Loud (RL) (Hooper et 
al. 1995) and Z97 showed that there are significant differ- 
ences in the strength of these lines between RL and RQ 
populations, independent of luminosity. A marginally sig- 
nificant RL/RQ difference was also seen in a sample of 255 
of the optically brightest QSOs in the LBQS sample (Fran- 
cis, Hooper, & Impey 1993). A more detailed examination 
of this effect will be made in a subsequent paper that will 
contrast the continuum and emission line properties within 
the full LBQS sample. 

In contrast, the Wx of the emission features in the F91 
composite are lower than the mean and medians of the 
measured emission lines in the LBQS sample. Although 
the F91 composite was created using a large proportion of 
the same spectra measured here, it is not surprising that 
the measurements should differ. Lyman a, which shows 
the largest dispersion of relative flux within the sample 
(sec Fig. 5 in F91), is most affected by narrow absorption 
features. Allowance for the effects of narrow absorption 
increases the Wx modeled with Gaussian components rel- 
ative to an estimate based on the integration of flux above 
a continuum. For [Ne III] and US, the F91 composite indi- 
cates a stronger feature than suggested by the KM mean 
presented in Table 7. However, the measurements are sim- 
ilar when only the detected lines are studied. 

There has been no investigation of inclusion of censored 
data on the emission line and continuum correlations, par- 
ticularly the Baldwin effect (Baldwin 1977; Zamorani et 
al. 1992; Osmer, Porter, & Green 1994) and the correla- 
tions associated with the primary eigenvector of spectral 
variance in AGN (BG92, M96, Sulentic et al. 2000). Pre- 
vious studies have focussed on a small number of bright 
emission features; particularly C IV A1549, He II A4686, 
H^, O III A5007 and iron emission. We briefly summarize 
the published measurements for QSO samples of BG92 (87 
low z QSOs, optical spectra), Mcintosh et al. (1999) (32 
high z QSOs with IR spectra), M96 (52 QSOs with optical 
and UV spectra) and Sulentic et al. 2000 (125 QSOs) in 
Table 8. Where upper limits are quoted, we have used sur- 
vival analysis to enable a more direct comparison to our 
results in Table 7. Some of the published work has also 
been presented separately as RL or RQ subsamples. While 
all these studies suffer from differing selection effects, the 
comparison samples also differ from the LBQS sample in 
redshift range and magnitude limits. 

The most notable difference between the properties of 
the LBQS and the measurements shown in Table 8 can be 
seen in the relative strength of optical iron emission, the 
mean value of which appears significantly higher in BG92 
and MS96. The region used to measure the Wx of optical 
iron emission is identical to that used in BG92 (and subse- 
quent authors) and a test using the methods described in 
§3 to measure the BG92 sample spectra (kindly supplied 
by T. Boroson) shows reasonable (^ 20%) agreement in 
Fe II Wx- The measurement of UV iron emission in M96 
appears similar to the mean of the distribution of iron 
emission strength measured in the LBQS sample. However 
the window used by M96 to characterize the iron emission 
strength was 1550 - I750A and the iron template flux in 
this region is approximately 6.5% of the flux in the region 
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used for the study of the LBQS sample presented here (i.e. 
the mean UV iron emission in M96 is 15 times higher than 
the mean in the LBQS sample!). The unusually strong 
iron emission was noted by M96 as a possible bias in their 
sample. The heterogeneous sample (52 AGN for which 
HST FOS spectra covering C IV were available in 1994 
and for which the authors had matching optical spectra) 
has an overabundance of strong iron emitters relative to 
the BG92 sample. The difference in iron emission distribu- 
tions contrasts to the similarity seen in the distributions of 
W\ for [O III]A5007. The strongest correlations reported 
for optical spectra of AGN involve iron emission, H/3 and 
[O III] . We expect that our new, more homogeneous emis- 
sion line measurements, of iron in particular, may signifi- 
cantly reshape the discussion of correlations that produce 
the largest variance in PCA analysis of quasar emission 
lines. 

7. SUMMARY 

This paper makes available measurements from the 
largest database of OUV spectra of QSOs to date. Our 
analysis of 993 spectra of similar quality from the LBQS 
yields over 8000 measurements of the 20 most prominent 
emission lines between rest frame 1025A - 5900A . We 



have accoimted for the effects of blended iron emission 
coniplexcis. narrow absorption lines, continuum breaks, 
and we include upper limits on the strength of low inten- 
sity emission lines. The measurements are presented with 
uncertainties generated in an objective manner and with a 
rigorous inclusion of censored data in tabulated emission 
line parameter distributions. A more exhaustive analysis 
of the relationships between the emission line and contin- 
uum properties of the LBQS sample will be presented in 
a later paper. 

High quality, reproducible measurements of emission 
lines for uniformly selected samples will advance our un- 
derstanding of active galaxies, especially in a new era of 
large surveys selected by a variety of complementary meth- 
ods, such as the Sloan Digital Sky Survey (SDSS; York ct 
al. 2000), the FIRST Bright Quasar Survey (White et al. 
2000), or the Chandra Multiwavelength Project (Green et 
al. 1999). 

The authors gratefully acknowledge support provided by 
NASA through grant NAG5-6410 (LTSA). CBF acknowl- 
edges the support of NSF grant AST 98-03072. KF thanks 
Matt Malkan and the UCLA Division of Astronomy & As- 
trophysics for their hospitality. 
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APPENDIX 

Notes on specific emission lines 

Lyman /? + O VI A1035 — A flat 'pseudo' continuum and a single Gaussian emission feature were used to model the 

spectrum near this emission line blend. The resulting equivalent widths should be viewed as approximate due to the 
nature of the QSO spectra in this region. Only in one spectrum of the sample was the Lyman /3 and O VI emission lines 
clearly distinct, but the spectrum was modeled with a single Gaussian for consistency. 

He II A1640 — The single Gaussian component used to model the region on the red side of C IV A1549 will account 
for omission from He II A1640, [No V] AAl.575,1593, [Nc IV] AA1602,1609, Si II AI650, [O III] AAI661, 1663, 1668, & Al II 
A1670. The use of multiple Gaussian components in all but a few spectra did not improve the model fit to this region due 
to the quality of the LBQS sample. The tabulated values for this emission line component, particularly the FWHM, will 
be much larger than expected for just the He II line. 

[Ne V] A3426, [O II] A3728, [Ne III] A3869 — The global continuum model was not successful for the region near 
these lines due to the blend of Balmer emission lines from transitions with m > 7 in the low resolution spectra and so a 
local power law continuum was created. This was measured from windows 30A wide at least 30A from the expected line 
position. The emission lines are measured above this local continuum. 

He II A4686 — For spectra where two Gaussian were used to model the H/J emission line, the H/3 broad component was 
included in the model fit of the He II A4686 emission line. 

Individual Objects 

LBQS 1206+1052 — This QSO was assigned a redshift of 0.402 ± 0.005 (Hewett et al. 1995) based on the cross 
correlation with the Francis et al. (1991) LBQS composite spectrum. Although all the features in the spectrum will 
contribute to the redshift estimate the stronger features will contribute more than the weak so it is unclear as to the 
cause of the discrepancy as the strongest lines of [O III] AA 4959,5007 are clearly not at the tabulated redshift. Only the 
measurement of the Mg II emission line agrees with this redshift, this is likely due to the presence of an absorption feature 
on the blue side of Mg II that is not modeled in the automated fitting procedure used here. Other lines in the spectrum 
all show a lower redshift. The value adopted here is ^; = 0.396 ± 0.003 based on the median of 10 emission measurements 
and the la value of the distribution of z. 

LBQS 0023+0228 — The spectrum of this QSO shows very strong and narrow (FWHM ~ instrumental resolution) 

forbidden emission lines, no evidence of broad H/3 emission and only a weak flat continuum. The strength of [O II] (W\ = 
148 til ^) and [O III] {Wx = 258 A) in the spectrum of LBQS 0023+0228 is much higher than aU the other QSOs 
in the LBQS except for LBQS 0004+0224 which shows a rising blue continuum and strong broad emission from Mg II 
and H/3 as well as having forbidden line FWHM > 1200 km s^^. The spectrum of LBQS 0023+0228 resembles that of a 
starburst galaxy (e.g. M82) rather than a QSO. The presence of these two strong narrow line objects in the sample does 
not aff'ect the parameter distributions due to the large numbers of objects in each sample. 

Electronic Tables 

The electronic version of Table 5 has an identical format to that published here but includes one row for each emission 
line that is listed in Table 3 even if the line lies outside the wavelength range of the observed spectrum. The electronic 
table contains zeros for unmeasured emission line parameters, rather than the (...) present in printed Table 5. This will 
aid in the use of this large machine readable ASCII table. The first line of the electronic table contains shortened column 
headings. This is followed by 31 lines for each QSO (a total of 30753 lines). The first line for each object tabulates the QSO 
designation and redshift and has a format (al0,f6.3). This is followed by 30 rows with a format (al6,6i7,3f 9.2,i3), 
note that the positive and negative error estimates appear in separate columns. Where only an upper limit for Wx of an 
expected emission line is measured, the value is given in the positive error column of Wx- The FWHM for these lines is 
given with zeros in the error columns (the position offset is by definition km ). Note that there are emission lines 
where the FWHM and position offset could not be constrained but the Wx is constrained to 2a. 

The format of the electronic version of Table 6 is identical to that presented here and as with the electronic version of 
Table 5, the positive and negative errors have separate columns. The first line of the electronic table contains the column 
headings and is followed by a single row for each object with the ASCH format (al0,2f 6.3,3f7.2,i6,3f 11 .5,3f7.2, iS). 
The normalization units for the power law continuum (columns 8, 9, and 10 in the electronic table) are 10^^** erg cm~^ 
g-i ^-1 stress that in the electronic form of Tables 5 and 6 zeros appear in all blank spaces that would normally 
appear blank. 
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Table 1 

lbqs qsos excluded from emission line measurements 



Designation 


Redshift 


r9i TV T 1 

Rer Notes 


Designation 


Redshift 


Ref 


Notes 


nnn4+ni47 


1.710 


n 


1 230+1 627R 


2.735 




1 2 


0010—0012 


2.154 


n 


1230+1705 


1.420 


III 




001 009Q 


2.083 


1 2 


1 9S1 +1 320 


2.380 


I 




0018—0220 


2.596 


1 


1235+1453 


2.699 


I 




nni 8+nn47 


1.835 


n 


1 235+0857 


2.898 


I 




nni Q+ni n? 


2.130 


n 


1 235+1 807R 


0.449 


I 




0021—0213 


2.293 


n 


1 239+0955 


2.013 


I 




0022+0150 


2.826 


n 


1240+1516 


2.297 




1 2 


0025—0151 


2.076 


n 


1240+1551 


0.573 


iii 


5 


002Q+001 7 


2.253 


n 


1240+1607 


2.360 


III 




0045—2606 


1.242 


V 


1242+1737 


1.863 




12 3 4 


0051— 001 Q 


1.713 


IV 


1243+0121 


2.796 


III 




0054+0200 


1.872 


IV 


1314+0116 


2.686 


III 




0059-2545 


1.955 


3,4 


1331-0108 


1.881 


III 




0059-0206 


1.321 


IV 


1332-0045 


0.672 




6 


0059-2735 


1.593 


V 


1333+0133 


1.577 


iii 


2 


0100-2809 


1.768 


V 


1442-0011 


2.226 


III 




0103-2753 


0.848 


V 


1443+0141 


2.451 


III 




0109-0128 


1.758 


IV 


2111-4335 


1.708 


V 




1009+0222 


1.349 


III 


2113-4345 


2.053 


V 


1,2 


1016-0248 


0.717 


III 


2116-4439 


1.480 


V 




1029-0125 


2.029 


III 


2140-4552 


1.688 


V 




1133+0214 


1.468 


III 


2201-1834 


1.814 


V 




1138-0126 


1.266 


III 


2208-1720 


1.210 


V 




1203+1530 


1.628 


III 


2211-1915 


1.952 


V 




1205+1436 


1.643 


I 


2212-1759 


2.217 


V 




1208+1535 


1.961 


I 


2241+0016 


1.394 


II 




1212+1445 


1.627 


I 


2244-0105 


2.030 


ii 


1,2,7 


1214+1753 


0.679 


I 


2350-0045A 


1.617 


II 




1216+1103 


1.620 


I 


2351+0120 


2.068 




1,2 


1219+1244 


1.309 


III 


2355-0209 


1.802 


II 




1224+1349 


1.838 


I 


2358+0216 


1.872 


IV 




1228+1216 


1.408 


I 











''Original LBQS paper that identified BAL (upper case) or associated narrow absorp- 
tion (lower case). 



Note. — 1) NAL on Lyman a; 2) NAL on C IV; 3) BAL near Lyman a; 4) BAL near 
C IV; 5) NAL on Mg II; 6) BAL near Mg II; 7) NAL on N V. 

References.— (I,i) Foltz et aL 1987; (II,ii) Foltz et aL 1989; (III,iii) Hewett et aL 
1991; (IV,iv) Chaffee et al. 1991; (V,v) Morris et al. 1991 
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Table 2 

Continuum and iron fitting windows 



Rest frame wavelength range (A) Emission lines nearby 

Continuum Iron Blueward Rodward 



1140 


- 1150"^ 






VI A1035 


Lyman a A 12 15 


1275 


- 1280 






N V A1243 


O I A1305 


1320 


- 1330 






O I A1305 


Si IV + O IV A1400 


1455 


- 1470 






Si IV + O IV A1400 


C IV A1549 


1690 


- 1700 






He II A1640 


Al HI A1859 


2160 


- 2180 


2020 


- 2120 


C III] A1909 


Mg II A2800 


2225 


- 2250 


2250 


- 2650 


3010 
3240 


- 3040*" 

- 3270 


2900 


- 3000 


Mg II A2800 


[Ne V] A3426 


3790 


- 3810 






[O II] A3728 


[Ne HI] A3869 


4210 


-4230 






m A4102 


H7 A4340 






4400 


- 4750*= 


[0 III] A4363 


H/3 A4861 


5080 


- 5100 


5150 


- 5500 


[0 III] A5007 


He I A5876 


5600 


- 5630 








5970 


- 6000 






He I A5876 


[N II] A6549 



''This window lies on the blue side of the Lyman a emission line and is only used 
where no other continuum window is available. 

''May have some iron emission contamination. 
°He II A4686 lies in this window. 
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Table 3 
Inventory of emission lines 





Ac 


FWHM'^ 


Window'' 




liiinissioii I1116 




yKiLL S ) 




Ly otes 


Lyman f3 A1025.7 + O VI A1035 


1030.0 


5000 


1010 - 1060 


See Appendix. 


Lyman a A1215.7 


1215.7 


7000s 


1170 - 1350 








2500n, 8000b 






N V A1241.5 


1241.5 


6000 






I A1305 


1305.0 


2500 






Si IV + IV A1400 


1400.0 


5800 


1350 - 1450 




C IV A1549 


1549.0 


6500s 


1350 - 1720 








3000n, 1 1000b 






He II A1640 


1640.0 


10000 




See Appendix. 


Al III A1859 


1859.0 


3500 


1820 - 1970 




C III] A1909 


1909.0 


5500 






Mg II A2800 


2800.0 


4000s 


2700 - 2900 


Iron emission may be strong in 






3000n, 8000b 




this window. 


[Ne V] A3426 


3426.0 


1000 


3390 - 3460 


See Appendix. 


[0 II] A3728 


3728.0 


600 


3700 - 3760 


Sec Appendix. 


[Ne III] A3869 


3869.0 


900 


3810 - 3930 


See Appendix. 


m A4101.7 


4101.7 


450 


4000 4200 




H7 A4340.5 + [O III] A4363 


4352.0 


3500 


4240 - 4440 




He II A4686.5 


4686.5 


1200 


4580 - 4790 


See Appendix. 


H/3 A4861.3 


4861.3 


4000s 


4750 - 5100 


Iron emission may be strong in 






lOOOn, 5500b 




this window. 


[0 III] A4959 


4959.0 


700 




Relative strength and separation of 


[O III] A5007 


5007.0 


600 




[0 III] lines is fixed in 1st iteration. 


He I A5875.6 


5875.6 


2000 


5825 - 5900 



^This is the FWHM used for estimating the upper limits of Wx in weak emission lines. More than one width is 
given for emission lines that can have two Gaussian components; s = single component, n = narrow component, b = 
broad component. 

''The wavelength range over which the emission lines are modeled. 
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Table 4 

Total number of emission lines modeled 



Fixed FWHM Wx TOTAL 

Emission Line and position'' upper limits out of 993 



UV iron 




294 


953 


Optical iron 




109 


247 


Lyman /3 + O VI 


10 


27 


130 


Lyman a single 


2 





164 


narrow 


U 


u 


yo 


broad 








96 


N V 


24 


8 


260 


I 


20 


121 


260 


Si IV + IV 


12 


19 


414 


C IV single 


3 


2 


408 


narrow 








80 


broad 








80 


He II A1640 


66 


69 


488 


Al III 


63 


181 


667 


cm] 


20 


26 


667 


Mg II single 


3 


42 


559 


narrow 








118 


broad 








118 


[Ne V] 


10 


365 


488 


[Oil] 


27 


272 


393 


[Ne III] 


12 


259 


363 


H 5 


6 


222 


309 


H 7 + [0 III] 


4 


87 


251 


He II A4686 


5 


155 


187 


H (3 single 


4 


18 


141 


narrow 








7 


broad 








7 


[0 HI] A4959 


6 


73 


148 


[O HI] A5007 


8 


54 


148 


He I 





30 


33 


TOTALS 


305 


2437 


8288 



''Does not include Wx upper limits. 
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Table 5 

Representative emission line measurements 



Designation (z) 












FWHM 


AAp 




Absorption 


Emission Line 


(km s"^) 


(km s"^) 


(A) 


Lines 


(1) 


(2) 


(3) 


(4) 


(5) 


0025+0009 z = 0.205 










UV iron 


000+9000 




701^^ 




Optical iron 


1200+^^;;'' 




7411^ 




Mg II single 


28001^«E;« 


-500ll« 






[Ne V] 


1200t^Jg° 


-lOOlgg 


c; r-l-lO.O 
^■^-4.3 




[0 II] 


210l^« 


100+140 
lUU_^4o 






[Ne III] 


1000+?r 


/^+2200 

'-'-2200 


5.01^4° 




H S 


5500lffog 


n+1100 

"-"-1000 


qcc-H5B 
oo_i4 




H 7 + [O III] 


2800+700 


''OU_3go 


O4-H20 
o"*_14 




He II A4686 


1200 




< 15 




H /3 single 


3000+^2^ 


300^240 


1051^^ 




[0 III] A4959 


6001280 




- - r\+40 




[0 III] A5007 


A rr\+120 

4501^^;; 




r,o+13 




He I 


6001!;^;; 




8-0_7.9 




2244+0020 z = 0.973 










UV iron 


4250^-^,° 




50.6li-^ 




Al HI 


64001^^0 


9001^1° 


12.8ltl 




CHI] 


5550_jgo 


-400l«°o 


27.9i:^-^ 




Mg II narrow 


40001^80 


-2501-r, 


35.6ltJ 




broad 


125001?^°" 


2001^°" 


27.4l«j 




[Ne V] 


1000 




< 2.8 




2354-0134 z = 2.211 










UV iron 


2000l«Z 




351^^ 




Lyman /3 + O VI 


540011™ 


isoolg" 


24.9l^:« 




Lyman a narrow 


3500lfo°o 


150i?°o 


33ilJ 




broad 


8100ltlo 






1 


N V 


54001^0 


n+160 

'-'-220 


34111 




01 


600iS° 


100+160 


1 3+1-5 


1 


Si IV + IV 


5000t™° 


n+380 
'-'-380 


13.4+3« 




C IV single 


eioolggo 


rO+180 


52.7+1? 




He II A1640 


17000+'^^'^'' 

l(UUU_3j^00 


-12001^^15 


26.0t^^g« 




Al HI 


3500 




< 7.0 




CHI] 


650011^0 


_40o+900 

WU_goo 


00+17 
"J^-12 




Note. — This is a digested form of the table available from the electronic jour- 



nal. Only emission lines measured in three example LBQS spectra are presented 
here, see §4 for more details. 



Emission Line Properties of the Large Bright Quasar Survey 



Table 6 

Representative continuum measurements 



Designation 
(1) 


Redshift 
(2) 


(3) 


(4) 


Power law continua 
Ac Ax 
(5) (6) 


Fa 

(7) 


Polynomial 
Order 
(8) 


0025+0009 
2244+0020 
2354-0134 


0.205 
0.973 
2.211 


3.010 
5.318 
3.368 


3.99 t\fi 

Z.OU _o.oi 
^•'i^J -0.24 


4220 
2170 
1463 


1 97 +0.16 
-0.15 

15.4 «;? 
1 iq +0.04 

J-.J-f -0.08 


1-08 tut 


2 
2 
1 



Note. — This is a digested form of the table available from the electronic journal. Only the 

continuum parameters from three example LBQS spectra arc presented here, see §4 for more 
details. The units for the normalization of the power law continuum {A\) presented here are 
10~^® ergs cm~^ s~^ A~^. In the electronic version of this table the normalization factor is 
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Table 7 

Emission line parameter distributions for the LBQS sample 



Wx FWHM 



Detected Kaplan-Meier Detected 



Emission Line 
(^'\ 


Num 


iviean 


iJLJ 


Num 


Limits 


Mean 
(7) 




Median 
\°) 


Num 

V') 


iviean 




Median 
K'-'-j 


UV iron 


659 


40.2 


23.0 


953 


294 


9Q Q 

zy.y 




0.8 


27.3 


659 


4610 ± 


130 


3470 


Optical iron 


138 


39.0 


21.5 


247 


109 


ZO.O 


_|_ 


1.6 


21.3 


138 


4630 ± 


310 


2600 


Lyman /3 + O VI 


103 


11.7 


10.1 


130 


27 


y .0 


HZ 


0.9 


7.3 


93 


4870 ± 


280 


4410 


Lyman a single^ 


164 


56.4 


28.8 


164 







-1- 


2.2 


53.0 


162 


7820 ± 


230 


7650 


narrow 


96 


28.1 


21.3 


96 




OS 1 
Zo.l 


zt 


2.2 


22.0 


96 


2650 ± 


120 


2470 


broad 


96 


68.3 


31.6 


96 




Oo.o 


_i_ 


3.2 


64.0 


96 


10250 ± 


420 


9330 


sum'' 


260 


71.2 


40.4 


260 





( L.I 


1 

± 


2.5 


61.0 










N V 


252 


18.8 


10.7 


260 


8 


1 o o 
lo.o 


1 

± 


0.7 


16.5 


228 


5580 ± 


170 


5400 


01 


139 


3.1 


2.5 


260 


121 


z.U 


_l_ 


0.1 


1.4 


119 


2990 ± 


190 


2460 


Si IV + IV 


395 


13.0 


8.8 


414 


19 


1 9 
IZ.U 


_l_ 
zlz 


0.4 


10.9 


383 


6780 ± 


160 


5940 


C IV single" 


406 


38.1 


19.3 


408 


2 


^s n 


_l_ 
zlz 


1.0 


34.5 


403 


7720 ± 


150 


7300 


narrow 


80 


17.8 


9.3 


80 




17.8 


± 


1.0 


15.5 


80 


2860 ± 


110 


2750 


broad 


80 


44.8 


23.1 


80 




44.8 


± 


2.6 


37.4 


80 


10960 ± 


360 


10530 


sum'' 


486 


42.1 


23.1 


488 


2 


42.0 


± 


1.0 


37.0 










He II A1640 


419 


18.7 


13.5 


488 


69 


16.4 


± 


0.6 


12.9 


353 


14670 ± 


390 


13590 


Al III 


486 


8.7 


6.0 


667 


181 


6.8 


± 


0.2 


5.2 


423 


5740 ± 


130 


5620 


C III] 


641 


28.4 


14.9 


667 


26 


27.6 


± 


0.6 


24.8 


621 


7820 ± 


170 


6670 


Mg II single" 


517 


39.1 


21.4 


559 


42 


37.1 


± 


0.9 


33.8 


514 


5160 ± 


120 


4440 


narrow 


118 


28.1 


12.9 


118 




28.1 


± 


1.2 


26.0 


118 


3510 ± 


110 


3370 


broad 


118 


36.6 


21.5 


118 




36.6 


± 


2.0 


31.6 


118 


8660 ± 


300 


8880 


sum*" 


635 


43.9 


25.1 


677 


42 


42.0 


± 


1.0 


36.8 










[Ne V] 


123 


5.8 


5.8 


488 


365 


2.1 


± 


0.2 


0.4 


113 


1570 ± 


90 


1420 


[Oil] 


121 


7.8 


14.6 


393 


272 


3.2 


± 


0.4 


1.0 


94 


900 ± 


70 


610 


[Ne III] 


104 


7.5 


5.0 


363 


259 


3.4 


± 


0.2 


1.8 


92 


1770 ± 


130 


1380 


m 


87 


12.3 


8.6 


309 


222 


3.8 


± 


0.4 


0.1 


81 


2860 ± 


220 


2300 


Uf + [O III] 


164 


18.3 


11.1 


251 


87 


14.6 


± 


0.7 


11.8 


160 


2920 ± 


140 


2500 


He II A4686 


32 


21.0 


20.1 


187 


155 


4.2 


± 


0.9 


0.1 


27 


3690 ± 


550 


2350 


H/? single" 


123 


62.4 


36.0 


141 


18 


57.4 


± 


3.1 


49.8 


119 


4370 ± 


250 


3760 


narrow 


7 


34.3 


13.6 


7 




34.3 


± 


4.8 


31.4 


7 


1160 ± 


270 


900 


broad 


7 


101.6 


58.9 


7 




101.6 


± 


0.6 


77.5 


7 


6560 ± 


850 


5300 


sum*" 


130 


66.4 


41.0 


148 


18 


61.2 


± 


3.4 


53.8 










[O III] A4959 


75 


13.8 


14.0 


146 


71 


8.5 


± 


1.0 


5.5 


69 


940 ± 


80 


770 


[0 III] A5007 


94 


30.4 


39.0 


146 


52 


21.6 


± 


2.8 


13.7 


86 


820 ± 


70 


590 


He I 


3 


10.7 


3.8 


33 


30 


4.7 


± 


1.1 


1.7 


3 


820 ± 


120 


750 



"The distribution for single Gaussian component models are tabulated separately from narrow and broad components. 

''The distribution of sum of the broad and narrow component Wx included with the single component measurements (see §5). 

Note. (1) Emission line or lino blend. (2) Number of detected emission lines, (3) Moan W\ of detected emission linos, (4) 
Standard deviation (SD) of W\ measurements for detected emission lines. (5) Total number of measured emission lines, (6) 
The number of upper limits estimates of Wx, (7)-(8) The KM reconstructed mean and median of the Wx distributions (see §5). 
(9)-(ll) The number, mean and median of the distribution of FWHM of the Gaussian components used to model each emission 
feature. All Wx are rest frame measurements in A and the FWHM are in km s~'. 
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Table 8 

Emission line parameter distributions for published comparison samples 



















FWHM 




Emission Line 


Rcf 


Num 


Limits 


Mcau 




Median 


Num 


Mean 


Median 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


UV iron 


I 


oo 
35 


o 
8 


oo 1 _1_ 

ZO.l ± 


A 

4.0 


1 K 

zl.o 








Optical iron 


1 

1 


0/ 


1 


4(.o ± 


o o 
2.0 


AO O 

43.0 










2 


45 


5 


178.5 ± 


6.2 


110.9 










3(4 


14 




25.6 ± 


4.9 












QT 
Olj 


1 1 




OA A _1_ 

34.4 ± 


A O 

4.Z 












4Q 


80 




52 ± 


2.8'' 












4L 


45 




21.5 ± 


2.5'' 










Lyman a 


3Q 
3L 


11 
6 




72.2 ± 
99.2 ± 


9.0 
9.0 










C IV single 


3Q 


13 




22.6 ± 


3.4 




15 


4880 ± 460 






3L 


13 




20.0 ± 


2.5 




15 


oOiO ± o70 




narrow 


o 
z 


o9 


U 


1 1 o _1_ 

ii.z ± 


1 A 

1.4 


5. 1 


on 
39 


1 ccn _l_ on 


1 ccn 


broad 


2 


52 





107.5 ± 


9.4 


96.0 


51 


6990 ± 310 


6610 


He II A1640 


2 


44 


3 


25.3 ± 


2.7 


23.4 


24 


12940 ± 600 


12400 


TTTl \ 1 rvrirv 
O iiij A19U9 




10 




17.9 ± 


1.6 




12 


6430 ± 730 






OT 


13 




1 ( .9 ± 


0. ( 




1 A 

14 


o/lon _l_ i?on 
04JU ± OzU 




TT \oonn 
Mg 11 AzoUU 


OT 


'7 
/ 

D 




Z3.1 ± 

on 1 _i_ 
Z9.1 ± 


o o 
Z.3 

A 1 

4.1 










Me 11 A4DOD 


1 
1 


I U 


u 


1 1 O 1 

11. z ± 


1.1 


i.O 










9 


oz 


o 


lU.O ^ 


z.^ 


lO.O 


14 


74sn + 77n 


7Qnn 


H/3 single 


1 


87 





95.5 ± 


4.0 


92.5 


87 


3790 ± 220 


3160 




3Q 


14 




64.8 ± 


6.9 




14 


4430 ± 590 






3L 


17 




73.2 ± 


5.1 




17 


5100 ± 470 






4Q 


80 




99 ± 


4.2'' 












4L 


45 




84 ± 


5.2'' 










narrow 


2 


42 





4.9 ± 


0.6 


4.0 


38 


630 ± 40 


590 


broad 


2 

3Q 
3L 


52 





105.3 ± 


8.4 


107.0 


51 
14 
17 


5950 ± 550 
9870 ± 970 
11890 ± 600 


4780 


[0 III] A5007 


1 


83 





24.8 ± 


2.7 


16.8 








2 


49 





58.8 ± 


6.8 


23.4 


49 


700 ± 40 


610 




3Q 


15 




20.6 ± 


3.5 




14 


1160 ± 90 






3L 


17 




10.2 ± 


1.9 




15 


1150 ± 130 






4Q 


80 




23.5 ± 


2.9'' 












4L 


45 




23.0 ± 


2.2'' 










This measurement 


is not suitable for direct 


comparison 


with the mean 


from the LBQS sample 


(Table 



8) because it is measured in a different region of the spectrum (see §6). 
''Errors on mean calculated from cr/\/iV. 

Note. — (1) Emission line or line blend. (2) Reference. (3) The total number of measured emission lines, 
(4) Number of upper limits on Wx. (5) (6) The mean and median of the W\ distributions. Where upper 
limits are present the KM estimator was used (see §5). (7)-(9) The number, mean and median of the 
distribution of published FWHM. All Wx are rest frame measurements in A and the FWHM are in km 

References. (1) Boroson & Green 1992; (2) Marziani et al. 1996; (3) Mcintosh et al. 1999; (4) 
Sulentic et al. (2000); Q = Radio Quiet, L = Radio Loud 
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Wavelength (A) 

Fig. 1.^ Continmim and iron template modeling windows. The composite QSO spcctrnm from Francis ct al. (1991) is 
shown above the iron emission templates for the UV (upper panel) and optical (lower panel) regions. The vertical dashed 
lines mark the positions of the narrow windows used to constrain the continuum shape for each spectrum. The spectral 
windows used in the modeling of iron emission are shown. See §§3.1, 3.2 and Table 2 for further details. 
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Observed wavelength (A) 

Fig. 2.— The continuum model for LBQS 1228+1116 (z=0.237). The bond between the polynomial continmim and the 
first power law is at 4000A in the observed frame and the inflection point between the two power law continua is at rest 
frame 4220A (see §3.1). 
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Fig. 3a.- The spectral models of three QSOs from the LBQS sample. For each QSO, the top panel shows 3 sections. The 
continuum model plotted over the observed spectrum including error bars on each bin, the continuum + iron emission 
template plotted over the spectrum, and the iron template profile alone. Note that for clarity the flux scales are different 
in these sections. Smaller panels for each emission line component show the total best-fit model plotted over the relevant 
region of each spectrum, the residuals, and the individual Gaussian components (and the profile of the iron template 
emission). Flux units are 10^^^ ergs cm^^s"^ wavelength units are in A and are observed frame values. Broad and 
narrow components of the same emission line species are marked 'b' or 'n' respectively. See §4 for more details, (a) — 
The spectral model of LBQS 0025+0009. 
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Fig. 3b.— The spectral model of LBQS 2244+0020 — See Fig.3a and §4. 
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Fig. 3c.— The spectral model of LBQS 2354-0134 — See Fig.3a and §4. 
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Fig. 4. The KM estimated distribution of tlic omission line properties of the LBQS sample. The npper panels are rest 
frame Wx (A) and the lower panels are the FWHM ( km s^^ ) of the Gaussian profile used to model the emission lines 
(see §3.2 and §4 for an explanation of the meaning of the iron emission Wx and FWHM). For emission lines that may 
be modeled with two Gaussian profiles, the narrow component distributions are marked (n) and broad components (b). 
The distribution of the sum of the narrow + broad component W\ is also shown as the dashed line in the panels for the 
single component models. 
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Fig. 4.— Cent. 
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Fig. 4.— Cont. 



